Quantum discord in spin-cluster materials 
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The total quantum correlation (discord) in Heisenberg dimers is expressed via the spin-spin corre- 
lation function, internal energy, specific heat or magnetic susceptibility. This allows one to indirectly 
measure the discord through neutron scattering, as well as calorimetric or magnetometric exper- 
iments. Using the available experimental data, we found the discord for a number of binuclear 
Heisenberg substances with both antiferro- and ferromagnetic interactions. For the dimerized anti- 
ferromagnet copper nitrate Cu(NC>3)2 • 2.5H2O, the three independent experimental methods named 
above lead to a discord of approximately 0.2-0.3 bit/dimer at a temperature of 4 K. We also deter- 
mined the temperature behavior of discord for hydrated and anhydrous copper acetates, as well as 
for the ferromagnetic binuclear copper acetate complex [Cu2L(OAc)] ■ 6H2O, where L is a ligand. 

PACS numbers: 03.65.Ud, 03.67.-a, 75.10.Jm, 75.50.Xx 



I. INTRODUCTION 

In the last years, it was understood that entanglement 
does not exhaust all quantum correlations in a system. 
Information correlations measured by entropy [lH3( were 
placed at the forefront. The total amount of correlations 
is identified with the mutual information I. The total 
correlations may be divided into two parts: purely clas- 
sical correlations C and purely quantum ones Q ■ 
The quantum excess of correlations, Q = I — C, has been 
called discord 

Note that the quantum discord has been evaluated 
explicitly only for several families of two-particle states 
(density matrices) 0-Q3- 

It is remarkable that discord can exist even in separa- 
ble (but mixed) states, i.e. when quantum entanglement 
is identically equal to zero. An example of two-qubit sep- 
arable state with nonzero discord Q — (3/4)log 2 (4/3) ~ 
0.311 was given in Ref. j 1 II ] - Below, in our paper, we 
will also discuss similar situations. As it turns out, "al- 
most all quantum states have nonclassical correlations" 
12]. Thus, quantum discord is a different measure of 
quantum correlation than entanglement. 

As shown on the model of deterministic quantum com- 
putation with one pure qubit (DQC1) [HI [Hj], quantum 
discord can lead to a speedup over classical computation 
even without containing much entanglement. Discord 
can also detect the quantum phase transitions [lj, [l5| . 
Moreover, recently it has been shown [l6| that in contrast 
to the entanglement and other thermodynamical quanti- 
ties, discord makes it possible to identify the quantum 
phase transition points at finite temperatures (T > 0). 

However, in order to utilize in practice the remark- 
able properties of quantum discord it is necessary to find 
methods to measure it experimentally. Our paper con- 
cerns this important problem. For the two-qubit Heisen- 
berg systems, we found the relations between discord and 



the ordinary spin-spin correlation function, as well as be- 
tween discord and internal energy, specific heat or mag- 
netic susceptibility. This allows one to determine the 
behavior of quantum discord in various substances and 
investigate its properties. 

The remaining part of this paper is organized as fol- 
lows. In Section 2, we describe the model and present 
formulas for calculating the total, classical, and quantum 
correlations. Here a comparison analysis is also made for 
those correlations and the spin-spin correlation function 
and for entanglement. Section 3 is devoted to expressions 
of discord in terms of spin-spin correlation function, in- 
ternal energy, specific heat, and magnetic susceptibility 
which are measured by standard experiments. These ap- 
proaches are used on the copper nitrate compound. In 
this section we also present the temperature dependen- 
cies for the discord in crystals of copper (II) acetate com- 
plexes, which are excellent examples of dimeric materials. 
The results obtained are briefly summarized in Section 4. 



II. CORRELATIONS IN HEISENBERG DIMER 

The Hamiltonian of a Heisenberg dimer reads 

1 



U 



(1) 



where J is the exchange coupling constant and = 
(af, cr?,af) the vector of Pauli matrices at the site i — 1 
or 2. The magnetic moment components for the dimer 
are 

M v = ~9»Vb(o-i +02), v = x,y,z. (2) 

Here, g v are the components of the g factor and /is is 
the Bohr magneton. 

The density matrix of a system in thermal equilibrium 
has the Gibbs form 
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P = ^cxp(-H/fc B T), 



(3) 
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where fcg is the Boltzmann constant and Z the parti- 
tion function which is found from the condition Trp = 1 . 
Performing necessary calculations for the system (JT]) one 
arrives at 



CORRELATIONS 



1 + G 

l'(T) = 7 (l+<»i*») = 7 I l 2G G 1-G 



with 



G(T) = -1 



3 + exp(-2J/k B T) 



1 + G 
(4) 



(5) 



It is easy to check that the quantity G equals the spin- 
spin correlation functions, 



G = (of of) = <o*o*> = Kaf), 



(6) 



where the brackets denote the statistical average. Values 
of G range from —1 to zero for the antiferromagnetic 
cluster ( J < 0) and from zero to 1/3 for the ferromagnetic 
one (J > 0). 

The density matrix (j3J) has the form for which the 
quantum discord is evaluated exactly Q : 

Q = I-C, (7) 

where the mutual information equals 

I = \ [(1 - 3G) log 2 (l - 3G) + 3(1 + G)log 2 (l + G)] (8) 
and the classical part of the total correlations is 

C = i[(l + |G|)log 2 (l + |G|) + (l-|G|)log 2 (l-|G|)]. (9) 
Expressions ([5]), ©-([5]) define so-called thermal discord 

[iMl- 

Another type of quantum correlation in a system is the 
entanglement of formation 
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where concurrence G for the ferromagnetic dimer com- 
pletely vanishes, and for the antiferromagnetic one is 

Bill 



G(T) 



-i(l+3G), T <T e 



0, T>T e . 
Here the temperature T e is given by equation 

1.8204.... 



\J\ e 



2 
ln3 



(11) 
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FIG. 1: The temperature dependencies of correlations \G\, 
Q, C, and E for the antiferromagnetic Heisenberg dimer. 



According to Eqs. (|71)- (fl"2")) , quantum discord and entan- 
glement are functionally related. 

In Fig. Q] the temperature dependencies of |G|, E, Q, 
and G, which describe different correlations in the anti- 
ferromagnetic dimer are shown. (In an antiferromagnet, 
the correlator G < therefore we took its absolute value.) 
It is seen that at T = all correlations are maximal (sat- 
urated)and equal to one. The correlations preserve prac- 
tically the same value with a small increase of tempera- 
ture. This is associated with the existence of a gap in the 
energy spectrum of the discussed system. With a further 
increase in temperature, all kinds of correlations mono- 
tonically decrease. By this, the curves pass through the 
inflection points where a convexity is changed to concav- 
ity. Functions \G(T)\, Q(T), and G(T) are different from 
zero for all temperatures T < oo. When T — > oo, the 
value of spin-spin correlations tends to zero as |G| ~ 1/T, 
and the information correlations Q{T) and C(T) go to 
zero more rapidly — according to the law 1/T 2 . 

From Fig.[T]it can be seen that at low temperatures the 
entanglement E is larger than both Q and G. Because 
E(T) has the sudden disappear point at the temperature 
T e , the curve E{T) must intersect both Q(T) and C(T). 
Coordinates of intersection points are 



k B T QE /\J\ 
Q(Tqe) - 



= 0.5880..., 
E(T QE ) = 0.7462 . . 



and 



k B Tce/ I J\ = 0.9260..., 
Q(T CE ) = E(T CE ) = 0.3390 . 



0.75 (13) 



(14) 



The phenomenon that in some situations the entangle- 
ment can be larger that the total quantum correlations 
has been pointed out in Ref. [7|. From our calculations, 
it follows that by T / discord is always greater than 
classical correlation, and the latter may be both larger 
and smaller than quantum entanglement. 

At the temperature T e , entanglement in the antiferro- 
magnetic dimer vanishes. On the other hand, the spin 
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FIG. 2: Correlations G, Q, and C vs ksT/J in the ferromag- 
netic Heisenberg dimer. 




FIG. 3: Discord as a function of G. Dashed vertical line 
separates antiferro- and ferromagnetic regions. 



correlation G(T e ) = —1/3, and the total entropy correla- 
tions 



I e = I(T e ) = 1 - i log 2 3 = 0.2075 



and discord 



Q(T e ) = 1 log 2 3 - ^ = 0.1258. 



(15) 



(16) 



Consider now a dimer with the ferromagnetic coupling 
J > 0. At zero temperature, G — 1/3 and the density 
matrix Q is 



Po = p(0) 



/ 2 

V 



1 1 
1 1 



= g[2(|00)(00| 



|11)(11|) + (|01) + |10))((01| + (10| 



(17) 



Here {|00), |01), |10), |11)} is the standard basis for the 
two-qubit system. The state (fT7|) is mixed because p§ ^ 
Po- At the same time, the state (|17p is separable (i.e., the 
entanglement E = 0). Indeed, a partial transposition of 
Po is 



p'd 



(18) 



in such a system. At zero temperature, the spin-spin cor- 
relation reaches its maximal value G = 1/3. Here quan- 
tum discord is also equal to 1/3. It is interesting that 
the total classical correlations are now less than Q and 
according to © equal only to 



Co = C(0) = - - log 2 3 = 0.0817 . 



(19) 
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Thus, the ratio of discord to classical correlation here 
achieves the value Q a /C = l/[5 - 31og 2 3] w 4.0798 (c.f. 
Ref. 0). 

At low temperatures all three types of correlations have 
quasi- horizontal sections ("pedestals"). With increas- 
ing temperature, the correlations pass through inflection 
points and then decay asymptotically to zero. Curves 
G(T), Q(T), and C(T) do not at any point intersect with 
one another. 

Notice the following. In accordance with Eqs. (0)-©, 
the quantum discord Q is a function of G (see Fig. [3]) . 
This function is monotonic both for the antiferromag- 
netic and ferromagnetic dimers. Moreover, if G = (ab- 
sence of ordinary correlations) then also Q — (absence 
of any quantum correlations), when \G\ is maximum (= 1 
or 1/3) then Q takes the same maximum values. This al- 
lows us to consider spin-spin correlation as a measure of 
discord. This is similar to that of the concurrence as a 
measure of entanglement [25J, [26J . 



Eigenvalues of this matrix equal 1,1,1, and 3. All 
these eigenvalues are positive. Consequently, in accor- 
dance with the positive partial transpose (PPT) crite- 
rion [iH, [23|], the state (flT|) is separable. Using formulas 
(jTj) — d9]> we find that the discord of a state (fT7|) equals 
Qo = 1/3 ~ 0.333. This value is larger than 0.311 for the 
discord of the example from Ref. [Il[ mentioned above. 

Let us look at Fig. [51 The figure shows the depen- 
dencies of G, Q, and C versus temperature for the dimer 
with the ferromagnetic coupling. Entanglement is absent 



III. EXPERIMENTAL MEASUREMENTS OF 
DISCORD 

As mentioned above, discord has a number of attrac- 
tive properties. Unfortunately, the question of how to 
measure it experimentally is open. But in the case of 
spin dimers, the information correlations /, C, and Q 
can be expressed via the experimentally observed char- 
acteristics of a system. 
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A. Neutron scattering and discord 

Inelastic scattering of thermal neutrons is a power- 
ful tool for the study of low-energy excitations in crys- 
talline transition-metal and other compounds. By this, 
the Fourier components of spin pair-correlation function 
are extracted from the scattering data f [27l - l29j and ref- 
erences therein). Performing the inverse Fourier trans- 
formation produces the correlation function itself. 

The neutron scattering experimental results for the 
quasi-dimer antiferromagnetic crystals of deuterated 
copper (II) nitrate CufNOg^ • 2.5D2O are presented in 
Ref. [30| (see also [3l|). (Notice that the deuteration 
needed for performing the neutron scattering experi- 
ments has no measurable effect on the exchange coupling 
32].) In the experiments with the copper nitrate, the 
neutron scattering intensity was measured for the tem- 
perature range 0.31 < T < 7.66 K. The temperature 
T e is equal to about 5 K [3l|. At T = 4 K, the spin- 
spin correlation has the value G = —0.54(9) [33j. Using 
Eqs. 0-© we estimate the discord as Q = 0.3(1). 

Note, in macroscopic systems, molar discord is Q = 
NaQ (Na is the Avogadro number). However, we will 
normalize discord per a dimcr. 



B. Internal energy and specific heat 

The internal energy per a mole of dimers ([1} is given 



u(T) = -^-G(T), 
2k H 



(20) 



where R = k B N^ is the universal gas constant. In turn, 
the energy equals 



u(T) =u + J c m (T)dT, 




(21) 



where c m (T) is the magnetic part of specific heat (the 
part after subtraction of lattice contribution from the 
total heat capacity). The integration constant uq in 
Eq. (f2~Tj) can be restored from the condition u(oo) = 0, 
i.e., 



u = tt(0) = - / c m (T)dT. 

" 



(22) 



Therefore, we can get the correlation function G, from 
calorimetric measurements, and then find the discord us- 
ing the expression of Q via G. 

On the other hand, the magnetic specific heat per mole 
of dimers is given as [3~i| 



c m (T) = 12R 



( J y exp(2J/k B T) 



c m /R 




FIG. 4: Relation of specific heat c m /R and correlation G. 
Dashed vertical line separates antiferro- and ferromagnetic 
regions. Points G = — 1 and G = 1/3 correspond to the 
temperature T — and G — to the T — 00. 



This function exhibits a maximum (a Schottky-like 
anomaly). For the ferromagnetic coupling, J > 0, its 
coordinates are 

k B T max /J = 0.9259..., CT/R = 0.1663..., (24) 

and for the antiferromagnetic one, J < 0, 

k B T max /\ J\ = 0.7029 . . . , C ax /E = 1.0234 .... 

Using Eq. (|5|) we can rewrite expression for the specific 
heat of Heisenberg dimers in the following forms: 

CmCn = ¥ (tV) (1 + G(T))(1 3G(T)) (26) 



4 \k B T r 
and, especially remarkably, 

c m /R = A (1 + G )(i _ 3G) In 2 (t^) ■ ( 27 ) 

These relations can be used to extract the correlation G 
directly from specific-heat measurements. The behavior 
of c m /R versus G is shown in Fig. |4] 

Consider as an example copper nitrate 
Cu(NC"3)2 • 2.5H 2 0. Its specific heat was measured 
for the temperature range 0.5 — 4.2 K [35| (see also [36[ 
and references therein). At temperatures below 0.5 K, 
specific heat is zero. Near 1.82 K it has the maximum 
value of 2 x 0.51i?. At higher temperatures, specific heat 
decreases obeying at T > 4 K the asymptotical law [35| 



l /2i? = 3.3/T 2 



(28) 



\k B T J 1 + 3exp(2 J/k B T) 



(23) 



The value of the exchange constant has been estimated 
as 2J/k B = -5.18 K |||. 

We have numerically integrated the available exper- 
imental data from zero to 4 K. Taking u(0)/R = 
3J/2k B = -3.885 K we determined that at T = 4 K 
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the internal energy is u(4)/R = —1.63 K and hence 
G(4) ~ -0.42 [see Eq. (J2DJ)]. Consequently, according 
to the calorimetric data, the discord in the given binu- 
clear cluster material at the temperature of 4 K equals 
Q = 0.19. This value can be considered to be reason- 
ably in agreement with the estimate obtained from the 
neutron scattering measurements. 

Note that the integration of function (j2"8")l from 4 K to 
infinity leads to u(A)/R = -1.65 K. Hence, G(4) ~ -0.42 
and we again return to the above result for the discord. 
On the other hand, according to Eq. (f2"8"| the specific 
heat at 4 K is c m (A)/R = 0.4125. Solving the transcen- 
dental equation (|27|) and taking into account that the 
temperature 4 K is larger than the maximum tempera- 
ture T max = 1.82 K, we find G ~ -0.4 (see Fig.|4j. Then 
Q ~ 0.18. 



C. Magnetic susceptibility 

Molar magnetic susceptibility of Heisenberg dimers 
satisfies the Blcaney-Bowers equation [3J, l37| 



2k B l 



(29) 



Here, g is the corresponding component of the Lande fac- 
tor when the measurements are made on a single crystal 
or 
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2 - o (5x + 9 2 v + 9l) 



(30) 



if the measurements are performed on a polycrystalline 
(powdered) sample. 

For the antiferromagnetic coupling (J < 0), the 
Bleaney-Bowers susceptibility displays a maximum with 
coordinates 



k B T a 



J 



1 + W(3/e) 



1.2472. 



[Jjw = l^(3/ e ) = 0.2011 



(31) 



(32) 



Here W(x) is the Lambert function defined by the 



equation We 



w _ 



This function under the name 



LambertW(x) was included in the Maple package. 
From Eq. ([2T?]) . we get the spin correlation function 



That is 



where 



G{T) _ ™> _ i. 



G(T) = -l + - X (T)/xo(n 



Xo(T) 



4,k B T 



(33) 



(34) 



(35) 



is the Curie law for the paramagnetic ions. Hence, one 
can indirectly measure the discord by performing mag- 
netometric measurements. 

Note that the Bleany-Bowers susceptibility is related 
to the internal energy ([20]) as 



u(T) = -3iV A J 



N A g^l X{T) 2 



Consequently, the magnetic specific heat is 



3k B J dT X (T) 



dT 



(36) 



(37) 



Moreover, inserting (|3"3")) into equation (I2T1) we obtain an 
expression for the specific heat c m /R trough the magnetic 
susceptibility without differentiation. 

We analyze, in the third instance, the sample with 
the copper nitrate compound. The magnetic susceptibil- 
ity of Cu(NC-3)2 ■ 2.5H 2 at low temperatures has been 
measured in Ref. (38|. The susceptibility of the powder 
has a rounded maximum at 3.2 K where it is equal to 
2 x 0.065 emu/mol. For the model of binary clusters and 
isotropic exchange, the authors [38( found from the data 
that -J/2k B = 1.28 K and g = 2.11. Using the exper- 
imental points presented in their figures (see Ref. [38() 
and taking into account our normalizations per a mole of 
dimers, we found that at T = 4 K the powder magnetic 
susceptibility is x — 2 X 0.063 = 0.126 emu/mol. Then, 
Eq. dH3j> yields G(4) = -0.396. As a result, the discord 
is Q = 0.17. This estimate agrees more or less with the 
values obtained above from the neutron scattering data 
and heat capacity measurements. Some discrepancy can 
be attributed to experimental errors and especially to the 
fact that the copper nitrate is only quasi-dimeric. 

We will now move on to other compounds, which be- 
long to the pronounced binuclcar materials. 



1. Antiferromagnetic substances 

Let us consider classical examples of spin-^ Heisen- 
berg dimeric materials crystals of copper (II) 
acetates, [Cu(CH 3 COO) 2 • H 2 0] 2 (hydrated) and 
[Cu(CH 3 COO) 2 ] 2 (anhydrous). Their magnetic suscep- 
tibility study has a long history beg inning in 1915 and 
continuing into the present day [39h43| . Experimental 
results between 90 and 400 K are given in table form in 
[40l | . The data have been described by the Bleany-Bowers 
equation with the fitted parameters 2J/k B — —408 K, 
g=2.13 for the copper(II) acetate monohydrate and 
2J/k B = —432 K, g=2.17 for the anhydrous copper(II) 
acetate [401 ] . 

Using Eqs. ([3"4"]) we obtained experimental 

points for the discord in both compounds. Results are 
presented in Fig.[S] Theoretical curves we plotted, taking 
the corresponding estimates for the coupling constants 
J/k B and using Eqs. ©, 
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Q, bit 




FIG. 5: Thermal quantum discord in the hydrated (-•-) and 
anhydrous ( — o — ) copper(II) acetates. 



of 1/3 



Q, bit 





FIG. 6: Quantum discord against the temperature in the 
ferromagnetic compound [Cu2L(OAc)] • 6H2O. 



We will now discuss the different domains in Fig. [5] 
According to Eq. (1131) the points which are higher than 
the level Q ~ 0.75 belong to the region where entangle- 
ment E > Q. However, at temperatures T > 120 — 127 K, 
the discord becomes larger than the entanglement in both 
compounds. The temperatures of 371 K and 393 K that 
are marked on the abscissa axis by the longer bars are the 
points of sudden disappearance of entanglement in hy- 
drated and anhydrous copper(II) acetates, respectively. 
At higher temperatures, the entanglement is zero. Near 
400 K, the discord is 11 — 12% of the maximum value 
Q(0) = 1. 



2. Ferromagnetic compound 

Entanglement in the ferromagnetic Heisenberg dimer 
is absent, but the discord is not zero. 

We found in the recent literature the exper- 
imental data for the ferromagnetic dimer mate- 
rial the binuclear copper (II) acetate complex 
[Cu 2 L(OAc)] • 6H 2 0, where H 3 L=2-(2-hydroxyphenyl)- 
l,3-bis[4-(2-hydroxyphenyl)-3-azabut-3-enyl]-l,3- imida- 
zolidine [44|. Magnetic susceptibility measurements for 
it have been performed between 5 and 300 K. Taking into 
account the crystal structure of this compound, the sus- 
ceptibility data have been fitted to the Bleany-Bowers 
equation. The best least-squares fit has been obtained 
with the parameters J/fcg 



35.4 K and g = 2.13 [44 1. 
At lowest measured temperatures near 5 K, the experi- 
mental points drop out from the theoretical dependence. 
This may be ascribed to the influence of weak interdimer 
couplings. 

Using Eqs. ©, and (J34J) we obtained from 

magnetic data the temperature behavior for the discord 
shown in Fig. [51 Maximal discord, Q ~ 0.32, is achieved 
near 13 K. At the temperature T = 300 K, the product 
\T = 0.89 cm 3 K/mol [44]. This leads to the discord 
Q ~ 0.003, i.e. about 1% relative to the theoretical limit 



IV. CONCLUSIONS 

We have related the quantum discord in binuclear spin 
clusters to the scattering data, as well as to the basic 
responses of a system to external perturbations — ther- 
mal and magnetic. This makes it possible to measure the 
discord in various magnetic substances. 

In this paper, we have done a comparison analysis for 
the temperature behavior of discord, classical correla- 
tions, entanglement, and spin-spin correlation function 
in the two-qubit Heisenberg systems. We have shown 
that in the case of antiferromagnetic interactions, entan- 
glement can be both larger and smaller than the discord 
or classical correlations. Discord and classical correla- 
tions can be present when entanglement is absent. This 
is observed both in antiferro- and ferromagnetic dimers. 
By T ^ 0, the discord is always larger than the classical 
correlations. 

We have presented the temperature dependencies of 
discord for different solid-state dimeric materials with 
antiferro- and ferromagnetic couplings. 

In the current literature, the transfer of quantum cor- 
relations from one system to another is intensively dis- 
cussed (e.g., p5H47| and references therein). In particu- 
lar, several schemes for extracting the entanglement from 
solids have been reported. With respect to discord, its 
dynamical behavior under the quantum state transfer, for 
example, between semiconductor double-dot molecules 
and photons is being studied now (48[. 
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